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In this paper we present a model for the short (< 1 second) population of 
gamma-ray bursts. In this model heated neutron stars in a close binary system 
near its last stable orbit emit a large amount of neutrinos (~ 10 53 ergs). A fraction 
of these neutrinos will annihilate to form an e + e~ pair plasma wind which will, 
in turn, expand and recombine to photons which make the gamma-ray burst. 
We study neutrino annihilation and show that a substantial fraction (~ 1/2) 
of energy deposited into e + e~ pairs comes from inter-star neutrinos, where each 
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Q_i' of baryon-free plasma that is deposited between the stars. To model the e + e 

pair plasma wind between stars, we do three-dimensional relativistic numerical 
hydrodynamic calculations. We find that the time scale for these bursts, deriving 
from the baryon-free plasma, is less than one second and they will have a hot 
spectrum ~ 5 MeV. The energy in bursts is the order of 10 52 ergs. 

Subject headings: gamma rays: bursts — gamma rays: theory 



member of the neutrino pair originates from each neutron star. Thus, in addition 
to the annihilation of neutrinos blowing off of a single star, there is a new source 



1. Introduction 

In Salmonson et al. (2001) we investigated a model for gamma-ray bursts (GRBs) deriv- 
ing from a neutrino burst of energy ~ 10 53 ergs above a heated, collapsing neutron star in a 
binary. Such conditions have been suggested by numerical relativistic hydrodynamic simula- 
tions (Mathews & Wilson 2000) of compression, heating and collapse of binary neutron stars 
near their last stable orbit. Such a thermal neutrino burst was found to partially recombine 
via uV — > e + e~ into an e + e~ pair plasma which expands relativistically. This fireball then 
recombines into photons via e + e~ — > 77 and was found to give a gamma-ray burst of energy 
~ 10 51 — 10 52 ergs with spectral and temporal characteristics consistent with observations. 
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This paper elaborates on the previous work in three distinct ways. 1) The latest general 
relativistic hydrodynamic calculations by Wilson and Mathews indicate that evolution and 
compression of the stars in a binary is faster than previously thought; < 1 sec, thus making 
this process a more natural candidate for the generation of "short" GRBs (< 1 sec). 2) 
Herein we use recent work by Salmonson & Wilson (2001) which calculates the neutrino 
annihilation rate between two neutron stars in addition to the annihilation rate from a single 
star (Salmonson & Wilson 1999). This 'inter-star' annihilation will not drive a baryon wind 
from the stellar surface and thus we find no baryon loading in the e + e~ pair plasma that 
originates between the neutron stars. 3) We employ three dimensional (3D) relativistic 
hydrodynamic calculations to model the flow of the wind in the complex, rotating, strong 
gravity environment around the neutron stars. 

It has been known for some time that the population of GRBs is bimodal, with ap- 
proximately a third of bursts having durations less than 2 seconds, and these short bursts 
typically have harder spectra than those bursts lasting longer than 2 seconds (Kouveliotou 
et al. 1993). This bimodality is thought to indicate that short and long bursts are separate 
populations and mechanisms. This view is bolstered by evidence that spectral break energy 
(Paciesas et al. 2001) and pulse lags (Norris et al. 2000) have a discontinous jump between 
short and long populations. The lack of discovery and location of an afterglow associated 
with short bursts has severely impeded progress in determining the distance scale and thus 
the energetics of these bursts. In fact, comparisons of BeppoSAX and BATSE data archives 
indicates a dearth of X-ray afterglows (Gandolfi 2001). Analysis of BATSE data for short 
burst decay tails give conflicting results indicating both a lack of an underlying afterglow 
component (Connaughton 2000) and the existence of such a component (Lazzati et al. 2001). 



2. The Model 

In this model we estimate that 10% of the rotationally enhanced binding energy of a 
neutron star (~ 10 53 ergs) is converted to thermal energy via compression, vortices and 
shocks. This 10 53 ergs of energy is released as a monotonically increasing luminosity of 
neutrinos over a timescale of order ~ 1/10 second as estimated by Mathews & Wilson 
(1997). This neutrino luminosity ~ 10 54 ergs s" 1 annihilates into an e + e~ pair plasma. 
Annihilation of high neutrino luminosities in the strong gravitational field of the neutron 
stars can have high efficiencies; near unity (Salmonson et al. 2001). About half of the pair 
plasma energy is deposited uniformly around the neutron stars due to single star neutrino 
annihilation (Salmonson & Wilson 1999) and the other half is deposited between the stars 
due to interstar neutrino annihilation (Salmonson & Wilson 2001). 
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This plasma deposition morphology then becomes input for the 3D relativistic hydro- 
dynamic code, which calculates the expansion of the plasma (Figure 1). These simulations 
show a plasma of very high entropy expanding out along the plane of symmetry between 
the neutron stars. In the regions around the stars lower entropy plasma is formed because 
a baryon wind is blown from the stars. The numerical model follows the baryons, hence 
some baryons were added in the interstellar region. The calculation yields an estimate of 
the angular distribution of the baryonic loading of the plasma. 

Thus this model predicts a variety of bursts. Viewed along the axis of rotation, a 
prompt quasi-thermal burst of duration ~ 1/10 second will result from the annihilation of 
fireball pairs (Salmonson et al. 2001). Because of the dearth of baryons left over to sweep 
into the interstellar medium, we do not predict the existence of an afterglow. This agrees 
with preliminary searches of the data archives for short-burst afterglows, which appear to 
be missing (Gandolfi 2001). 

Viewed far from the axis of rotation a very different burst results. The lower entropy 
means that there will not be a prompt burst from pair annihilation in the fireball. However, 
there will be a baryon wind sweeping into the interstellar medium. Thus we expect a burst 
that decays into an afterglow as a power-law (Salmonson et al. 2001). This behavior will be 
made chaotic and complex by the rapid rotation of the binary system. 

3. One-Dimensional Calculations 

Using one-dimensional (ID) relativistic hydrodynamics described in Salmonson et al. 
(2001), we can study the expansion of an e + e~ pair plasma fireball from the surface of a 
single neutron star. As per our model, we deposit £ = 10 52 ergs of energy into a narrow 
(width = 1 km) shell region above the surface of a Rq = 10 km radius and 1.4 M mass 
neutron star. We monotonically ramp the energy over At = 0.1 second. This model is run 
with a range of deposited baryon densities and thus a range of baryon loadings of the fireball. 

It is important to note that our hydrodynamics do not track the protons and neutrons 
separately and, as such, assumes they are well coupled during the expansion and acceleration 
of the plasma. As was shown by Fuller et al. (2000) and Derishev et al. (1999), the neutrons 
can decouple from the protons and e + e~ plasma during acceleration, thus changing the 
overall baryon load of the plasma. Since it is plausible that the baryon matter blown off of 
the neutron star is 90% neutrons, it is reasonable to expect, depending on the efficiency of 
the decoupling, the effective baryon mass to be a factor of several less than the total baryon 
mass. 
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The plasma fireball is evolved and we find that it expands adiabatically as a thin shell 
of constant coordinate thickness AR (Ruffini et al. 2000). We track this expansion and pair 
annihilation until the optical depth, r, of the plasma falls below unity: r ~ Na^AR < 1 
where N is the coordinate density of electrons and positrons and ot is the Thompson cross- 
section. The energy in photons, ^photons , and in the baryons, -Ebaryons, is calculated and these 
results are summarized in Fig. (2). One can see that for £/M < 10 5 , where M is the baryonic 
rest mass energy, most of the energy goes into accelerating the baryons while for £/M > 10 5 
a substantial fraction of the energy is expressed as photons directly from the plasma; i.e. a 
gamma-ray burst. 

This result can be understood when one considers that the opacity of the fireball derives 
from two electron densities: electrons associated with entrained protons, N e i ectrons = N protons , 
and electrons and positrons in pairs, N pairs . Defining the optical depth as r = r e i ectrons + 
Tpairs ~ (N pro t ons + N pairs )(j T AR, then by comparing the densities N protons and N pairs , we 
understand the relative importance of these two sources of opacity. 

We first consider the baryons. Equating the energy deposition with outgoing flux gives 
the initial energy density at the surface of the neutron star E = E/ (AIA^Rqc) ^3x 10 29 
ergs cm -3 . It follows that the corresponding initial number density of baryons is iV^o = 
Eq l 'rripC 2 (M / '£) ~ 1.3 x 10 35 (M/£)cm -3 , where M/£ is the ratio of baryon rest mass energy to 
total energy. Finally, we know (Ruffini et al. 2000) that for relativistic adiabatic expansion we 
have: iVj, = Nb t0 (R /R) 2 . Let us assume that 10% of the baryons are protons with associated 
electrons. Then the optical depth from the entrained baryons is T e i ectrons ~ Q-IN^UtAR ~ 
0.1(E /m p c 2 )(M/ 8,)(Ro/R) 2 a T AR and thus optical thinness, T e i ectr0 ns ~ 1 happens at a 
radius R ba ryons ~ 10 12 y/(M/8.)- 5 cm, where (M/£)_ 5 = (M/£)/10- 5 . Notice that this 
radius depends on the initial baryon loading of the fireball: Rbaryons oc a/M/£. 

The opacity due to pairs is a bit more complex. As discussed in Salmonson et al. (2001), 
the pair number density is governed by the following equation 

= -~C-^N pairs Vn + a^((N; mrs (T')) 2 - N 2 mrs )/W 2 (1) 

where o ann v is the Maxwellian averaged mean pair annihilation rate per particle and N pairs (T') 
is the local coordinate equilibrium e + e~ pair density at local temperature T' given by the 
appropriate Fermi integral with a chemical potential of zero. Early in the fireball ex- 
pansion e + e~ pairs are in equilibrium so N pairs N pairs (T') and thus the pairs expand 
adiabatically: N pairs = N pairS £(Ro/R) 2 , where the initial pair number density is given 
by N pa i r . : o ~ (ai?Q) 1 / 4 /3A; ~ 10 34 cm -3 . When the local temperature, T', of the plasma 
falls below T[ ~ m e c 2 /3 ~ 1/6 MeV the pairs begin to rapidly annihilate. This hap- 
pens at a radius R* ~ R /T'(E / 'a) 1 / 4 ~ 100.Ro- The pair number density plummets 
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until the annihilation term in Eqn. (1) no longer dominates, which happens at a density 
Npair,* ~ R*/ (cr ann Rl) ~ 10 20 cm~ 3 . This annihilation happens over a few which is very 
small compared to the radius at which pairs and photons decouple and therefore these re- 
sults are independent of the details of this annihilation phase. Beyond this point the pair 
number density again expands adiabatically: N pa i r = N pa i r ^(R*/R) 2 . Thus the optical depth 
from the pairs is T pairs w N pairs a T AR ~ {E /a) 1/A a T AR/ (og^RoT'l)(R /R) 2 and optical 
thinness, T pairs ~ 1, happens at a radius R pa i rs ~ a/ (E /a/Tl) 3 / A a T RQAR/o\^ ~ 10 12 cm 
which is independent of £/M. 

Thus we see that the transition shown in Fig. (2) from baryon-dominated to photon- 
dominated fireball energy is governed by the relative importance of opacity associated with 
the baryons and that associated with the e + e~ pairs. Therefore it is important to consider 
both the opacity due to electrons associated with protons, which dominate for high baryon 
loading, and the opacity of the e + e~ pair plasma which dominates for low baryon loads. 
For initial values £/M < 10 5 , -Rbaryons > -Rpairs so the opacity is dominated by the electrons 
associated with the baryons and thus the baryons are able to "hold" the plasma longer and 
thus absorb all of the fireball energy. For £/M > 10 5 we have i? pa irs > -Rbaryons and the pairs 
are the dominant source of opacity and thus determine the radius of optical thinness. If the 
radius of optical thinness falls below the radius to saturate the energy of the baryons, i.e. the 
radius at which most of the plasma energy is in the baryons, -R S aturate = (£/M)i?o, then the 
acceleration of the baryons is inefficient, they do not attain the energy (£/M)m p c 2 , and thus 
more energy remains in the photons, thus yielding a photon burst. 



While the ID hydrodynamics of the last section elucidate the behavior of an expanding 
relativistic plasma, a 3D calculation is required to capture the structure inherent in the 
present problem of rotating, emitting binary neutron stars. In this section we introduce the 
3D relativistic hydrodynamic equations. Herein we use units G — c — 1. The continuity 
equation is 



where the spatial dimensions are indicated by the sum over i = 1, 2, 3. The energy equation 
is 



4. 



Three-Dimensional Calculations 
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Fig. 1. — Three dimensional relativistic simulation of two 10 km radius neutron star sepa- 
rated by 30 km emitting 10 53 ergs/sec of energy in e + e~ pair plasma and about 1 % equivalent 
mass in baryons. The contour map, with right star cutaway, is of baryon density. The vector 
field is the 3-velocity of expanding plasma. This problem settles down to a static flow after 
about one orbit with period ~ 1/300 second. 
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Fig. 2. — Output of the ID code run over a range of initial energy-to-mass ratios, £/M, 
with energy fixed at £ = 10 52 ergs. This shows that the energy available for creation of an 
internal burst versus that for an external burst varies as a function of baryon loading. The 
more baryon loading (lower £/M), the more energy is entrained in the baryons, -Ebaryons, 
and thus creates a more energetic external shock. Less baryon loading (higher £/M) makes 
a more efficient fireball burst, where much of the original energy is left in the photons, 
-E'photons- The middle figure shows the typical observed temperature, T Q b s , of the resultant 
photons. The lower figure shows the characteristic radii determining the evolution: -Rbaryons ~ 
10 12 a/ (M/ £)_5 cm, -Rp a i r s ~ 10 12 cm, R sa turate ~ 10 6 (£/M) cm. The photons decouple from 
the baryons at the maximum, -Rbaryons or -Rpairs- If Rsaturate is greater than either Rbaryons or 
Rpairs, then the baryons are unable to saturate with energy and thus the energy remains in 
the photons. 
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And the momentum equations are 



dS„ = _ 6S dJo^ _ _ JP _ W(D + ft. 



dt 



dt 6 dx* 



dx i 



dl3 j 



W - 



1 

w 



<9(log < 



where the four- velocity is 



(4) 



Ui 



D + YE ' 



and we define a generalized Lorentz factor 



W = W1 + 



A ' 



the coordinate velocity in terms of the four velocity (eqn. 5) is given by 

alii 



V 1 = 



W(p 4 



and the equation of state is given by 



r = 1 + 



PW 



(5) 



(6) 



(7) 



(8) 



The metric is 3+1 and taken to be conformally flat. The lapse function is given by 



a 



1 - ^ 

2_ 

1 + 1 



(9) 



The spatial factor is 



and the shift vector is 



= 1 + 



2 ' 



(10) 



oc Q x r — A I 



Mi M 2 
ri r 2 



(11) 



The factor F is defined by 



F = 



M _|_ M 

n r2 



Ml 3 r{ 



ri > ^ 



(12) 



where r\ and r 2 are the distances from the center of the nearest and farthest stars respectively, 



each of mass M and radius R*. 
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Fig. 3. — The terminal (final) Lorentz factor, 7fi na i = 2'jE/D, for ^^EjD ^> 1, from 3D 
simulations as a function of angle from the rotation axis of the binary system in the plane 
described by said axis and the line connecting the two stars. The local Lorentz factor is 7. 
Thus 0° is along the rotation axis and 90° corresponds to being directly behind one of the 
stars. The most energetic material, 7fi na i — > 00, is located within a narrow region, 9 < 15°, 
on the symmetry plane between the two stars. 
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Fig. 4. — From the 3D code, the observed temperature at infinity, T obs oc a(W 3 E) 1 ^ of the 
baryon-free pair plasma along the z-axis of rotation. Most of the plasma is at the high value 
of z hence this observed average effective temperature will be near 5 MeV. 
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5. Results 

An example of a 3D run is shown in Fig. 1. An edit of the final Lorentz factor, 7fi na i, as 
a function of angle from the rotation axis in the plane defined by the rotation axis and the 
line connecting the two stars is shown in Fig. 3. One can see that the highest energy per 
baryon occur along the rotation axis. From this figure one gets the scale of the opening angle 
of this high energy material: ~ 15°. The simulation does not sufficiently resolve the region 
of high E/D at low angles. Very little baryon mass will make it into this region of high 
energy density, thus we expect E/D to be effectively infinite on the rotation axis. Thus for 
viewers within this angle we have a thermal burst dominated by e + e~ pair-photon plasma as 
described in (Salmonson et al. 2001). An edit of the observed temperature T t, s oc a(W 3 E) 1 ^ 4: 
is shown in Fig. 4, where we have used quantities defined in eqns. (3, 6, 9). 

One can see from Fig. 3 that, depending on the observer angle with respect to the axis 
of rotation, a wide range of baryon loadings, parameterized by E/D, will be observed. Thus 
we study the range of expected observed signals as it depends on E/D. Herein we assume 
a correspondence between the ratio of energy to mass, £/M, of the spherically symmetric 
calculations of Section 3 to the ratio of energy to mass coordinate densities, E/D, as taken 
from the 3D code runs. As such, we expect a thermal burst for observers along the axis of 
rotation of the stars and within an angle ~ 15° where no baryons are deposited and thus 
E/D — > oo. As seen in Fig. 3, for increasing 9 from the rotation axis, E/D decreases rapidly 
to ~ 100. 

For all E/D = £/M < 10 5 all of the energy of the pair plasma is deposited into the 
baryons and therefore a prompt thermal burst does not occur. Instead, the mass, M, of 
baryons moving with Lorentz factor, 7 final = £/M, sweep into the interstellar medium, thus 
forming a collisionless shock, and decelerate according to (e.g. Piran 2000) 

£ = 7finaiM = 47T P 7 2 r 3 (13) 



for 

/ o \ 1/3 

vL6 -2/3 F l/3 -1/3 



47r P c2 7final 



r>r dec =[- ~ 10 lb 7 ^ a E& V ' cm (14) 



where r dec is the deceleration radius, p is the interstellar medium (ISM) density, n\ = p/m p 
cm" 3 is the ISM baryon number density, 7300 = 7/300 and £52 = £/(10 52 ergs). This 
model was worked out in detail in Salmonson et al. (2001). The shock will emit synchrotron 
radiation as it decelerates with a characteristic luminosity L ~ (n^l^^-K^lO 51 ergs s" 1 
at a characteristic photon energy e ~ 200y / ni7| 00 keV and a duration on the deceleration 
timescale t dec = r dec /2^ 2 c ~ 3 (_E 52 /n 1 /7| 00 ) 1 / 3 s. Note that for 7 < 300 this timescale is 
consistent with long bursts: > 3 s. 
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A distinct feature of this progenitor model, since the wind comes from the surface of a 
neutron star, is that most of the baryons will be in the form of neutrons (say 90%). Thus 
these neutrons will not couple with the ISM magnetic field and will pass by the decelerating 
protons. These neutrons will decay back to protons on a timescale of a neutron half-life, 
r n,i/2 ~ 1000 s, which happens at a radius Rdecay = 7cr ni i/ 2 ~ 7300IO 16 cm. Comparing 
this radius with the decelaration radius, r^ ec , we see that r dec > Rdecay for 7 < 300, i.e. the 
neutrons have decayed into protons before reaching the deceleration radius. So we expect 
this effect to be important for 7 = £/M > 300 (Fig. 5). 

Another characteristic Lorentz factor 7 = £/M is when the time, in the fluid frame, to 
reach the deceleration radius is shorter then the deposition time, At = 0.1 s. The fluid frame 
deposition distance is Ar' = ^cAt ~ 7 300 1 12 cm. Also in the fluid frame, the deceleration 
distance is r' dec = r dec /^ ~ 10 15 (ni7f 00 )~ 1//3 cm - These are equal at 7 ^ 1500 (Fig. 5). Thus 
for 7 > 1500 the plasma deposited at the beginning of the deposition time has begun to 
decelerate before the plasma has completed deposition, so the assumption that the plasma 
is a thin, shocked shell plowing into the ISM is invalid. Instead, energy will be continuously 
deposited into the shock. 

An important caveat for these high Lorentz factors, 7 3> 300, is that they likely are 
beyond the regime of validity for the collisionless shock. The detailed physics behind such 
shocks is not well understood (e.g. Piran 2000) and assumes that the kinetic energy of the 
baryons, (7 — l)m p c 2 , will effectively amplify the magnetic field and energize the electrons. 
However, low densities of ultra-high energy baryons will inefficently couple with the ISM and 
thus the collisionless shock approximation will be invalid. 

Again, for energy, £ = 10 52 ergs, with photon energy e ~ 2007| 00 keV, the observed 
photon fluence seen at a distance, R, will be ££/47nR 2 e ~ 250£7^o^Gp C cm_2 ) where Rc pc = 
R/{1 Gigaparsec) and £ is an efficiency factor which may be about 10%. For example, given 
7 = 1500, we expect a fluence of 0.04 photons cm -2 with characteristic energy e ~ 0.1 GeV 
over a duration t dec 0.04 s. So such a burst would be far from recognizable as such (Fig. 5). 



6. Discussion 

A key result of the 3D numerical simulations (Fig. 3) is that the emission from this 
system is bimodal: about half of the total energy is deposited as a pure, baryon-free, E/D — > 
00, pair plasma along a 'fan' of angular half-width 6 = 15° along the symmetry plane between 
the neutron stars, and the other half of the total energy is blown off of the neutron stars as 
a wind with E/D m 300. Very little of the energy is deposited in intermediate regimes of 
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Fig. 5. — A schematic figure representing the physical regions along the spectrum of 7 = £/M 
discussed in Section 5. A simple external shock results for 7 < 300. For 300 < 7 < 1500 the 
neutrons will be decoupled from the protons and thus will decay and shock at a larger radius, 
thus creating a smeared external shock. For 1500 < 7 < 10 5 the external shock cannot be 
considered a simple thin shock and the theory of collisionless shocks is likely violated. In the 
region 7 > 10 5 , the energy remains in the photons from the original pair plasma and thus a 
simple internal shock results as discussed in Section 3 and is detailed in Fig. (2). 
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E/D. As such we analyse the expected observations from this model. 

The solid angle subtended by the fan of baryon-free plasma is fi/ an = 47rsin6 l and the 
wind blows off into the rest of space, Vl wind = An (I — sin 9). Therefore the ffuence from each 
component at a distance R for total energy £ is 



££_ 

(15) 



>n n fan R 2 8nsin9R 2 



7 _ ££/2 _ ££ 

wind 



& wind R 2 8vr(l - sin#)i? 2 

where £ is an efficiency factor. Now the system is rotating so an observer with angle, if), from 
the rotation axis will see a mixture of the fan and wind, unless he is located within an angle 
ip < 9. Since the period of rotation of the system is much smaller than the energy deposition 
timescale, At = 0.1 s, we average over the fluences from the wind and the fan to get a total 
average fluence. The proportion of a rotation subtended by the fan is 2/7rarcsin(#/ sin^) 
for 9 < ip ^ tt/2 and is unity for < ip < 9. So the average fluence contributed by the fan 
as a function of viewing angle is (Fig. 6) 

8%R I rib for0<^<# , 



and that of the wind is 



1— 2/-7T arcsin(0/ sin(i/>)) 
1— sin t 



for < ^ tt/2 , 



8vri? 2 [ for < ip < 9 . 

The total average fluence is 

Fty) = F /an (^) + F^ nd (^) . (18) 

An observer within < ip < 9 sees only the fan emission. This is effectively a jet of opening 
half-angle 9 = 15°. Such an observer sees a jet of strong thermal emission over a timescale 
of 0.1 s. For 9 = 15° = vr/12 we have 



F(rl> = 0) = 1.6 x 10~ 4 ergs cm" 2 

Rr 



F(ip = it/2) = 7.3 x 10~ 5 ergs cm 

^Gpc 



£-^5 2 
i? 2 



(19) 



so an observer along the jet (ip < 9) will see roughly twice the fluence in a short, 0.1 s, 
thermal, T obs ps 5 MeV, gamma-ray burst, than an off-axis observer (ip ~ 7r/2) who will see 
a synchrotron burst of duration ~ 3 s with characteristic photon energy 200 keV. 
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Fig. 6. — The top figure shows the average fluence F of the fan (Eqn. 16) and the wind 
(Eqn. 17) as a function of viewing angle ip. The bottom figure shows number fluence, F 
divided by the characteristic photon energy for the fan, Zf an ~ 5 MeV (Fig. 4), and the 
wind, e wind 200 keV (Sec. 5). This demonstrates that most of the energy comes from the 
fan while most of the photons come from the wind. A viewer within ip < 15° = n/12 will 
see only the fan emission, which effectively constitutes a jet of half-opening angle 15°. 
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In conclusion, when the binary system is viewed within 15° of the axis of rotation, we 
have a model for short gamma-ray bursts, ~ 0.1 s, which have hard spectra, several MeV. 
In addition, when the binary system is viewed at angles larger than 15° from the rotation 
axis, our model yields a second group of bursts, 30 times more common, with soft spectra, a 
few hundred keV, and intermediate time duration of 1 to 5 s. The energy available depends 
on the neutron star equation of state and the masses of the neutron stars. The energy in 
a burst could range from 10 50 ergs to a few times 10 52 ergs. From Fig. 6 we see that if the 
detector is primarily sensitive to number counts then at intermediate viewing angles the fan 
contribution to the signal may be missed. The pure fan signal, i.e. a jet viewed within 15° 
of the rotation axis, will have no afterglow, but the fan plus wind would have an afterglow 
and a small contribution of high energy photons. Jensen et al. (2001) have observed a 2 s 
burst with an afterglow. Unfortunately the high energy detector data was not obtained. 

This work was performed under the auspices of the U.S. Department of Energy by 
University of California Lawrence Livermore National Laboratory under contract W-7405- 
ENG-48. 
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